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Abstract 

In this paper, some initial-boundary- value problems for the time- 
fractional diffusion equation are first considered in open bounded n- 
dimensional domains. In particular, the maximum principle well- 
known for the PDEs of elliptic and parabolic types is extended for the 
time-fractional diffusion equation. In its turn, the maximum principle 
is used to show uniqueness of solution to the initial-boundary-value 
problems for the time-fractional diffusion equation. The generalized 
solution in sense of Vladimirov is then constructed in form of a Fourier 
series with respect to the eigenfunctions of a certain Sturm-Liouville 
eigenvalue problem. For the one-dimensional time-fractional diffusion 
equation 

(D?u)(t) = — (p(x)—)~q(x)u + F(x,t), a:G (0,1), t € (0,T) 



dx \ dx 

the generalized solution to the initial-boundary-value problem with 
the Dirichlet boundary conditions is shown to be a solution in the 
classical sense. Properties of the solution are investigated including 
its smoothness and asymptotics for some special cases of the source 
function. 
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1 Introduction 

During the last few decades, PDEs of fractional order have been successfully 
used for modeling of many processes and systems that stimulated develop- 
ment of their mathematical theory. This process is still going on, but one 
can already recognize some new insights and correlations between different 



branches of mathematics and natural sciences that have emerged within the 
framework of FC models. The literature related both to theory of PDEs of 
fractional order and to their applications is very extensive. For the theory 
of fractional differential equations we refer e.g. to the books [2], [9], and 
|24j . As to the applications, we mention the book [25] that is completely 
devoted to different applications of FC in physics, chemistry, technique, as- 
trophysics, etc. and contains several dozens of interesting case studies. The 
recent book [3] deals with applications of fractional differential equations 
in classical mechanics, group theory, quantum mechanics, nuclear physics, 
hadron spectroscopy, and quantum field theory. For other interesting FC 
models we refer the reader to [3], [5], [7], [TO], [12], [TO], [TO], [20], [22] to 
mention only few of many recent publications. 

In applications, relevant processes like e.g. anomalous diffusion or wave 
propagation in complex systems mainly run within some bounded domains 
in space that corresponds to initial-boundary-value problems for the PDEs 
of fractional order, which model the processes under consideration. Both 
analytical and numerical investigation of these problems thus become an 
important task in Fractional Calculus. 

In the author's papers |11|-|15|. the time-fractional diffusion equation 

(Dfu)(t) = div(p(x) grad u) — q{x)u + F(x,t), < a < 1 

with the Caputo fractional derivative Df with respect to the time variable 
t along with some of its important generalizations (multi-term equation and 
equation of distributed order) and some initial and boundary conditions has 
been considered in open bounded n-dimensional domains. In particular, 
the maximum principle well-known for the PDEs of elliptic and parabolic 
types was extended for the time-fractional diffusion equation as well as for 
the multi-term time-fractional diffusion equation and for the time-fractional 
diffusion equation of distributed order. The validity of the maximum prin- 
ciple is based on an important extremum principle for the Caputo fractional 
derivative. In its turn, the maximum principle was used to show uniqueness 
of solution to some initial-boundary-value problems for the time-fractional 
diffusion equations. This solution - if it exists - continuously depends on the 
problem data, i.e., on the source function and on the initial and boundary 
conditions. 

As to the existence of solution, in the author's papers [TT], [TO], [TO], the 
Fourier method of variables separation was used to obtain a formal solution 
in form of a Fourier series with respect to the eigenfunctions of a certain 
Sturm-Liouville eigenvalue problem. Under certain conditions, the formal 



solution was shown to be a generalized solution in sense of Vladimir ov, i.e., a 
continuous function that satisfies the PDE of fractional order in generalized 
sense (see [22]). To prove that this generalized solution is a solution in the 
classical sense, i.e., it is at least twice differentiable function with respect to 
the spatial variables and a-differentiable with respect to the time variable 
turned out to be a difficult undertaking that is still not resolved in general 
case. 

In this paper, this problem is solved for the initial-boundary-value prob- 
lems with the Dirichlet boundary conditions for the one-dimensional time- 
fractional diffusion equation. Under certain conditions, the unique solution 
to this problem can be represented via a Fourier series with respect to the 
eigenfunctions of a corresponding Sturm-Liouville eigenvalue problem. 

For other methods and results in the theory of initial-boundary-value 
problems for the PDEs of fractional order we refer e.g. to [I], [6], [8], |21j . 
[23] . |25| to mention only few of many recent publications. 

The rest of the paper is organized as follows. In the 2nd section, def- 
initions and notions used in the paper along with the problem formula- 
tion are given. In the 3rd section, main results for the initial-boundary- 
value problems for the time-fractional diffusion equation in open bounded 
n-dimensional domains are short presented. In particular, the maximum 
principle well-known for the PDEs of elliptic and parabolic types is extended 
for the time-fractional diffusion equation and used to show uniqueness of 
solution to the initial-boundary-value problems for the time-fractional dif- 
fusion equation. The generalized solution in sense of Vladimirov is then 
constructed in form of a Fourier series with respect to the eigenfunctions of 
a certain Sturm-Liouville eigenvalue problem. The 4rd section deals with 
the one-dimensional time-fractional diffusion equation. In this case, the gen- 
eralized solution is shown to be a solution in the classical sense. Properties 
of the solution are investigated including its smoothness and asymptotics 
for some special cases of the source function F. Finally, the last section 
contains some conclusions and open problems. 

2 Problem formulation 

In this paper, we deal with the time-fractional diffusion equation 

{D«u){t) = -L{u) + F{x,t), (1) 



< a < 1, (x,t) e Q T := G x (0, T), G C R n , 



and especially with its one-dimensional case (n = 1), where 

L(u) := — div(p(x) grad u) + q(x)u, 

p G C^G), g G C(G), p(s) > 0, q(x) > 0, s G G, (2) 
the fractional derivative Df is defined in the Caputo sense 

(B a f)(t):=(I 1 - a f')(t), 0<a<l, (3) 
I a being the Riemann-Liouville fractional integral 

'rfeft-^/W^ o<a<i, 

/(t), a = 0, 



and the domain G with the boundary S is open and bounded in R n . For 
the theory of the fractional integrals and derivatives the reader is referred 
e.g. to [9j or [21]. 

For a = 1, the equation ((TJ) is reduced to a linear second-order parabolic 
PDE. The theory of this equation is well-known, so that the main focus 
in the paper is on the case < a < 1. For the equation ©, the initial- 
boundary-value problem 

u\ t=0 = u (x), x G G, (4) 
u| s = v(x,t), (x,t) G S x [0,T] (5) 

is considered. A solution of the problem flTJ), ©, © is called a function 
it = u(x,t) defined in the domain fix ■= G x [0, T] that belongs to the 
space C(n T ) n Ty/((0,T]) n C|(G) and satisfies both the equation (HJ) and 
the initial and boundary conditions (jlj)-©. By W/((0, T]), the space of the 
functions / G G 1 ((0,T]) such that /' G L((0,T)) is denoted. If the problem 
CO)) ®) © possesses a solution, then the functions F, uq and v given in 
the problem have to belong to the spaces C(0,t), C(G) and C(S x [0,T]), 
respectively. In the further discussions, these inclusions are always supposed 
to be valid. 



3 n-dimensional time-fractional diffusion equation 

In this section, a short survey of results for initial-boundary-value prob- 
lems ©-© for the n-dimensional time- fractional diffusion equation CO) i s 
presented. For the proofs of results we refer the readers to |13j and |14j . 

We start with uniqueness of solution of the problem © , @ , © . The 
main component of the uniqueness proof is an appropriate maximum prin- 
ciple for the equation ©. 



Theorem 1 Let a function u G C(Qt) n W/((0,T]) n C%(G) be a solution 
of the generalized time-fractional diffusion equation (T7P in the domain Qt 
and F(x, t) < 0, (x, t) G Q T . 

Then either u(x,t) < 0, (x,t) G &t or the function u attains its positive 
maximum on the part S£ := (G x {0}) U (S X [0, T\) of the boundary of the 
domain Qt, ie-, 

u{x, t) < max{0, max u(x, t)\, V(x, t) G Oy. (6) 
(x,*)esg 

Similarly to the case of the PDEs of parabolic type (a = 1), an appro- 
priate minimum principle is valid, too. 

The maximum and minimum principles can be applied to show that the 
problem ([I]), ([H)-© possesses at most one solution and this solution - if it 
exists - continuously depends on the data given in the problem. 

Theorem 2 The initial-boundary-value-problem ^-$$) for the time- 
fractional diffusion equation (QP possesses at most one solution. This so- 
lution continuously depends on the data given in the problem in the sense 
that if 

\\ F ~ ^1lc(fi T ) - e > ll n - Uo\\c(G) - e 0> \\ v ~ v \\c(Sx[0,T]) < e l> 

and u and u are the solutions of the problem (QP, (^-(EP with the source 
functions F and F , the initial conditions uq and u~q, and the boundary con- 
ditions v and v, respectively, then the norm estimate 

r-pa 

h - u\\c(Q T ) < max{e , ei} + - ^ e (7) 
for the solutions u and u holds true. 

To tackle the problem of existence of the solution, the notion of the 
generalized solution is introduced (see [29] for the case a = 1). 

Definition 1 Let F k G C(fir), u ok G C(G) and v k G C(S x [0,T]), k = 
1, 2, . . . be the sequences of functions that satisfy the following conditions: 
1) there exist the functions F, uq, and v, such that 

\\ F k ~ F\\ c( p T ) ^0 as +oo, (8) 

\\uok ~ u o\\c(G) ~> as k +oo, (9) 



IK - v\\ C (Sx[0,T\) ~+ as k ->• +00, (10) 

#) /or any fe = 1, 2, . . . i/iere exists a solution u k of the initial-boundary-value 
problem 

Uk\ t=0 = uok(x), x eG, (11) 
u k \ s = v k {x,t), (x,t) eSx [0,T], (12) 
for the time-fractional diffusion equation 

(D?u k ){t) = -L{u k )+F k (x,t). (13) 

Suppose, there exists a function u G C(Or) such that 

\\ u k — u \\c{G) ~ * as ^ ~~ ^ +°°- (14) 

TTie function u is called a generalized solution of the problem |7J), 0)-(0). 

A generalized solution of the problem (pQ), ©-([5]) is a continuous func- 
tion, not a generalized one. Still, the generalized solution is not required to 
be from the functional space C(Q:r)nW/((0, T])nC%(G), where the solution 
has to belong to. 

It follows from Definition 1 that if the problem ([TJ), (J2])-(j5]) possesses a 
solution then this solution is a generalized solution of the problem, too. In 
this sense, Definition 1 extends the notion of solution of the problem ([T]), 
([H)-©. This extension is needed to get some existence results. But of course 
one does not want to lose uniqueness of the solution. The following result 
is valid: 

Theorem 3 The problem |7]), (2P-(EP possesses at most one generalized 
solution in the sense of Definition 1. The generalized solution - if it exists 
- continuously depends on the data given in the problem in the sense of the 
estimate |7j). 

Contrary to the situation with solution of the problem (pQ), (p[|)-(j5|), ex- 
istence of the generalized solution can be shown in the general case under 
some standard restrictions on the problem data and the boundary S of the 
domain G. In this section, existence of the generalized solution of the prob- 
lem 

(Dfu)(t) = -L(u), (15) 
u \t=o = u o( x )' x £ G, (16) 
u\ s = 0, (x,t) <e S x [Q,T] (17) 



is considered to demonstrate the technique that can be used with the ap- 
propriate standard modifications in the general case, too (see the case n = 1 
in the next section). The generalized solution of the problem (|15p -(|17p can 
be constructed in an analytical form by using the Fourier method of the 
variables separation. Let us look for a particular solution u of the equation 
(|15p in the form 

u(x, t) = T(t) X(x), (x, t) G Qt, (18) 

that satisfies the boundary condition f)lT[) . Substitution of the function (|18p 
into the equation (|15p and separation of the variables lead to the equation 

(D?T)(t) L(X) 

~ ~Mx) ~ ' ( } 

A being a constant not depending on the variables t and x. The last equation, 
together with the boundary condition (|17p . is equivalent to the fractional 
differential equation 

(D^T)(t) + \T(t) = (20) 
and the eigenvalue problem 

L(X) = XX, (21) 

X\ s = 0, x G S (22) 

for the operator L. Due to the condition (|2|), the operator L is a positive def- 
inite and self-adjoint linear operator. The theory of the eigenvalue problems 
for such operators is well-known (see e.g. [29]). In particular, the eigen- 
value problem (|2ip - (122p has a counted number of the positive eigenvalues 
< Ai < A2 < • • • with the finite multiplicity and - if the boundary S of G 
is a smooth surface - any function / € Ml can be represented through its 
Fourier series in the form 

00 

f{x) = Y J U,X, i )X i (x), (23) 
i=i 

where AQ € Ml are the eigenfunctions corresponding to the eigenvalues A,: 

L(Xi) = X t X u i = 1,2,.... (24) 



By Ml, the space of the functions / that satisfy the boundary condition 
(p2j) and the inclusions / G C 1 ^) n C 2 (G), L(f) G L 2 (G) is denoted. 



The solution of the fractional differential equation (|20p with A = Aj , i = 
1,2, .. . has the form (see e.g. [16], [17] ) 

Ti(t) = Ci E a {-U a ), (25) 

i? a being the Mittag-Leffler function denned by 

^W-E^FFij- <26) 

fc=l v ; 

Any of the functions 

u i (x,t) = c i E a (-X i t a )X i {x), i = l,2,... (27) 
and thus the finite sums 



(x,t) = J2c i E a (-X i t a )X i (x), k = 1,2... (28) 



i=l 



satisfy both the equation (|15p and the boundary condition (|17p . To construct 
a function that satisfies the initial condition (|16p . the notion of a formal 
solution is introduced. 

Definition 2 A formal solution of the problem U5\) -i71) is called the 
Fourier series in the form 

oo 

u[x, t) = Y^ua,Xi)E a (- Xit a ) Xi(x), (29) 

8=1 

Xi, i = 1,2,... being the eigenf unctions corresponding to the eigenvalues Aj 
of the eigenvalue problem h21\) - and (/, g) the standard scalar product 
in L 2 (0,l). 



Under certain conditions, the formal solution (J29J) can be proved to be 
the generalized solution of the problem f)15[) - (j 1 T[) . 

Theorem 4 Let the function uq in the initial condition i!6\) be from the 
space Ml. Then the formal solution I129\) of the problem &15\)-[T7\) is its 
generalized solution. 



To prove that this generalized solution is a solution in the classical sense, 
i.e., it is at least twice differentiable function with respect to the spatial 
variables and a-differentiable with respect to the time variable turned out 
to be a difficult undertaking that is still not resolved in general case. 

In the next section, this problem is solved for the initial-boundary-value 
problem for the one-dimensional time-fractional diffusion equation. 

Remark 1 The method presented in this section can be applied with some 
small modifications to the case of the problem M5\ )-(T7\) on an infinite do- 
main !] = G x (0, oo), G C R n , too. 

4 One-dimensional time-fractional diffusion equa- 
tion 

All results presented in the previous section are valid for the n-dimensional 
time- fractional diffusion equation and thus for the one-dimensional case, too. 
The main result of this section is that the generalized solution of the initial- 
boundary-value problems for the one-dimensional time-fractional diffusion 
equation we received in the previous section can be interpreted as solution 
in the classical sense under some suitable conditions. 

We start with the initial-boundary-value problem for the inhomogeneous 
one-dimensional time-fractional diffusion equation 



(Dfu)(t) = — [p(x)—)-q(x)u + F(x,t), x G (0, 1), t G (0, T) (30) 




with the inhomogeneous Dirichlet boundary conditions 



u (o,t) = Mt), 



u(l,t) = cj> 2 (t), 0<t<T, 



(31) 



along with the initial condition 




(32) 



and the conjugate condition 



n (0) =<M0). 



(33) 



In the following, we always suppose that the conjugate condition (|33|) holds 
true. Like in the case of the initial-boundary-value problems for the PDEs of 
the parabolic type (a = 1 in (J3QJ)), the problem (|30"|) - (f33|) can be transformed 



to a problem with the homogeneous boundary conditions. Indeed, let us 
introduce an auxiliary function v according to the formula 

v{x,t) = u{x,t) + j{Mt) - <h(t)) - Mt)- (34) 



Then the inhomogeneous boundary conditions (|3ip for the unknown function 
u are transformed to the homogeneous boundary conditions for the function 
v: 

v(0,t) = 0, v(l,t) = 0, < t < T. 



The form of the time-fractional diffusion equation (|30p for the new unknown 
function v 

d ( dv \ 

(D?v)(t) = — (p(x)—)-q(x)v + F 1 (x,t), x e (0,1), t € (0,T) 
with 

F 1 (x,t) = F(x,t) + j(D?(<Pi ~ fa))(t) - - Mt) ~ Mt) p'(x)+ 

q(x) (j(Mt) - Mt)) - Mt) 

as well as the form of the initial condition (1321) 



v(x,0) = u (x) + j(Mo) - Mo)) - MO) 



remain the same as for the original problem f|30|) - (j32|) . 

Without loss of generality, we can thus consider the initial-boundary- 
value problem p0j) - ([32|) with the homogeneous boundary conditions, i.e. 
with 

Mt) = Mt) = o, 0<t<T. 

Let us start with the homogeneous equation ([30]) . i.e. we first suppose that 
F(x, t) = 0, < x < I, < t < T. In the previous section, it was shown (see 
Theorem 4) that if uq 6 M l then the formal solution (|29p is the generalized 
solution to the problem (|30l) - (|32l) . We prove now that (f29j) is a solution in 
the classical sense under some additional conditions. 

Theorem 5 Let uq € A4l and L(uq) 6 Ml- Then solution of the initial- 
boundary-value problem 



lt=o 



u Q {x), < x < I, (35) 



u(0,t) = u(l,t) = 0, < t < T 
for the one- dimensional time-fractional diffusion equation 



(36) 



{D>){t) = ^ { P[X) ^) ~ q[x)U (3?) 

exists, is unique, and is given by the formula 

oo 

u(x,t) = ^2(u ,Xi)E a (—Xit a ) Xi(x), (38) 
i=i 

Xi, i = 1,2,... being the eigenf unctions corresponding to the eigenvalues Xi 
of the eigenvalue problem 121\) - Ii22\) . 

Proof: The uniqueness of solution follows from Theorem 2. Theorem 4 
states that (|38|) is a generalized solution to the problem (f35]) - ([3"T|) . We show 
now that (|38|) is a solution in the classical sense, i.e. it is at least twice dif- 
ferentiable function with respect to the spatial variable and a-differentiable 
with respect to the time variable. Indeed, let us differentiate the series from 
the left-hand side of the formula (|38p with respect to the spatial variable 
term by term and construct a series of derivatives: 

oo 

J^dEai-^nxKx), c i = (u ,X i ). (39) 
1=1 

Because of the inclusion uo G Ml, the function L{uq) is from the space 
L 2 (0, 1). We can thus determine the Fourier coefficients of L(uq) with respect 
to the eigenfunctions Xi, i = 1,2,...: 

(L(u Q ),Xi) = (uo,L(Xi)) = Xi{u Q ,Xi) = X iCi . (40) 

In the formula (j40f) . we used the fact that L is a positive definite and self- 
adjoint operator. The Parseval equality for L(uq) leads to the relation 

oo 

]TA^ 2 = ||L(uo)|| 2 <+oo. (41) 

i=l 

It is well-known that the system Xi, i = 1, 2, . . . of the eigenfunctions of the 
eigenvalue problem (I2ip - (122p is a complete orthonormal system in L 2 (0, 1) 
and the series 



are uniformly convergent on the interval [0,1] (see e.g. [29]). 

Using this fact, the estimate (|4ip . and applying the Cauchy-Bunjakovski 
inequality we can now prove the uniform convergence of the series ([39]): 



Y J \ciE a {-\ i t a ) X[{x)\ < MY,\ci\ 
i 

2/00 



i=l 



8=1 



< 



M E« E 



vi=l 



,i=l 



A? 



1 

|2\ 2 



< 



M\\L(u Q )\\ [J2 



,i=l 



A? 



|2 \ 2 



In the first inequality, we used the estimate (see e.g. [24]) 

M 

\E a (-x)\ < < M, 0<x, < a < 1. 

1 -)- x ~ 



(43) 



Because the series ([39]) is uniformly convergent one, too (see [13J ) , we thus 
have proved that the generalized solution ([38]) is a Cj(0, 1) function and the 
relation 

' " ^2*E a (-\i1?)Xl(x) 



dx 



i=l 



holds true. The same arguments can be applied to show that ([38]) is a 
Cj(0, 1) function and the relation 



d 2 u 
dx 2 



J2^E a (-\it a )Xl\a 



(44) 



8=1 



holds true. Indeed, due to the inclusion L{uq) € Ml we get L 2 (uq) € £ 2 (0, £) 
and 

(L 2 (u ),Xi) = (u ,L 2 (Xi)) = \ 2 (u ,Xi) = \ 2 a, i = 1,2, ... . 



Then 

J2 x t c i = \\ l2 ( u o)\\ 2 < +00 
1=1 

and we can use the uniform convergence of the 3rd of the series in (|42|) and 
again apply the Cauchy-Bunjakovski inequality to prove the formula 



i=l 



i=l 



XV (x) 



3 

A? 



< 



Note, that in the 2nd inequality we used the inequality Af < CXf, i = 
1,2,... with a certain constant C that holds true because of lim^oo \ = 
+00. 

To prove that the series given by (|38|) is a-differentiable with respect to 
the time variable we use the following result from |26j : 
Let for a sequence of functions /j, z = 1,2,... defined on an interval (a, b] 
following conditions be fulfilled: 

1) For a given a > there exists fractional derivatives (D a fi)(t), t £ (a, 6] 
for all i = 1, 2, . . . , 

2) both the series Xa=a/*(^) ana - ^ e ser i es X)i=iC^ a /i)(*) are uniformly 
convergent on the interval [a + e, 6] for any e > 0. 

Then the function defined by the series * s a-differentiable and the 

relation 

00 00 

p a £/i)(t)=£(i?7i)(t), a<t<6 
i=l i=l 

holds true. 

Now let us apply the fractional derivative of order a with respect to the 
time variable to the series from the left-hand side of the formula (|38|l term 
by term and construct a series of derivatives: 

00 00 

J2c i (D a E a (-\ i t a ))(t)X i (x) = -^c J A^ Q (-A 4 t Q )X i (x). (45) 

i=l i=l 

The formula ()45|) follows from the fact that the Mittag-Leffler function 
E a {—\it a ) is a solution of the fractional differential equation ()20p with A = 
Aj, ? = 1,2,... or can be directly verified using the well-known formula 

(0 < a < 1, < p) 

< D ' Vfl >< f > " fiT 1 ^^" - < 46 » 

Then we apply the Cauchy-Bunjakovski inequality to the series (|4"5j) and get 
the estimates 



J2\ -c i X i E a {-X l t a )X l {x)\ < M^|qA 
i=i i=i 



X l {x) 
1 

A, 2 



< 




MWEefAH [J2^IL\ < Ml \\L\u )\\ ± lX ^ 



i 

2 \ 2 



\i — l I — \ ' / \ I — J ' / 

Once again, in the 2nd inequality we used the inequality Xf < CXf, i = 
1,2, . . . with a certain constant C. It follows now from the uniform con- 
vergence of the 1st of the series in (|42p that the series (|45p convergences 
uniformly on any interval [e,l], < e < /. The function (|39p is thus a- 
differentiable and the formula 

oo 

(D?u)(t) = -Ys^iM-U^Xiix) (47) 

i=i 

holds true. 

We have shown that the function defined by the series (|38p is at least 
twice differentiable function with respect to the spatial variable and a- 
differentiable with respect to the time variable and thus is the solution in 
the classical sense. ■ 

Let us now discuss smoothness of solution (I38p with respect to the time 
variable. The following result is valid: 

Theorem 6 Let uq € Ml and L(uq) E M-l- Then the solution ( fggj) of 
the initial-boundary-value problem l^35\) - (3b}) for the one- dimensional time- 
fractional diffusion equation j37| ) is a continuously differentiable function 
with respect to the time variable on the interval (0,T). 

Proof: We already know (see Theorem 5) that (|38p is the solution of the 
problem (j35l) - (|36l) for the equation (|37l) in the classical sense. Let us differ- 
entiate the series from the left-hand side of the formula (|38p with respect to 
the time variable term by term and construct a series of derivatives: 

oo 

Y t -CiX i t a ~ 1 E a>a (-X i t a )X i (x), a = (u ,X i ), (48) 

i=l 

E a> p being the generalized Mittag-Leffler function defined by 

OO 

**>«-EfjdE+?r (49) 



To get the expression (|48p . the formula 
d 



^_E a (—Xit a ) — —Xit a 1 E a ^ a (— Xit a ) 



was used. The last formula can be obtained by the term by term differen- 
tiation of the power series for the Mittag-Lemer function and by using the 
well-known properties of the Gamma-function. 

Now let t be from the interval [e,T], < e < T. Using the estimate (see 
e.g. [2U) 

M 

\E a ,p{-x | < — — , 0<x, 0<a< 1, 0< p, (50) 
1 + x 

we get the inequality 

M \t a M 
lA^- 1 ^ J-Xit a )\ < — < — < C, e < t < T (51) 

with a constant C that of course depends on e but not on t. 

Then we use the estimate (|5ip and the same reasoning as in the proof of 
Theorem 5 and get the following chain of estimates: 



i=i i=i 



XAx] 



i 



< 




<h sftiiK«.)ii E 



\Xi(x)f 



i i \ i • •' / \ i X{ J 

vl=l / \l = l / \« = 1 / 

In the 2nd inequality, we used the inequality Aj < C2Xf , z = 1,2,... with a 
certain constant C%. It follows now from the uniform convergence of the 1st 
of the series in (j42[) that the series (|48|) uniformly converges on any interval 
[e,T], < e < T. The function (f39|) is thus continuously differentiable for 
< i < T and the formula 

^-^CiAjrXot-Vl^W- Ci = («o,^) (52) 

i=l 

holds true. ■ 

Now we consider the case of the initial-boundary- value problem (|35p -(j36l) 
for the inhomogeneous one-dimensional time-fractional diffusion equation 
([30]). 

In the following, we suppose that the initial condition uq and the function 
L{uq) are from the space Ml- In addition, we suppose that the source 
function F belongs to the space Ml for any t € (0, T), too. This means, 
that for any t € (0, T) the function F can be represented in form of an 
uniformly convergent Fourier series with respect to the spatial variable x: 



oo „i 
F(x,t) =Y / F i (t)X i (x), Fi(t) = (F,Xi) = / F(x,t)Xi 
i=i J o 



x) dx, (53) 



Xi, i = 1, 2, ... being the eigenfunctions corresponding to the eigenvalues Xi 
of the eigenvalue problem ([2"T]) - (f22j) . 

As in the case of the PDEs of parabolic type (a = 1 in (|30p ). we look 
for solution of the initial-boundary- value problem ([35]) - ([36]) for the equation 
(|30p in form of a Fourier series 

oo 

u{x,t)=Y J T i {t)X i {x), T i (t) = (u,X i ). (54) 
i=i 

From (|54|) , the Fourier series for the initial condition (remember the inclusion 
u G Ml) 

oo 

uq(x) = s ^c i X i {x), Ci = (u ,Xi), 

i=l 

and the uniqueness of the Fourier series for a function from the space Ml, 
we get the initial values for the unknown Fourier coefficients 

T l (0) = c i = (u ,X l ), z = 1,2, (55) 



Substituting now the representation (|54j) into the equation (|30j) and us- 
ing the properties of the eigenfunctions Xi, we get an uncoupled system 
of initial- value problems for the fractional differential equations for the un- 
known Fourier coefficients Tj, z = 1,2,...: 

{(D«T i ){t) = -\T l {t) + F l (t), 
\Tm = Ci. 

The solution of the initial- value problem (|56|) is known (see e.g. [16] . 
[17J) and is given by 

Ti(t) = Ci E a (-\t a ) + f r a ~^ E a ^ a (—XiT a )Fi(t - t) dr. (57) 

J o 

Thus the formal solution to the initial-boundary- value problem (|35|) - (|36|) for 
the inhomogeneous one-dimensional time-fractional diffusion equation ()30|) 
can be written in the form 



u[x 



t)=jt (ciE a (-\it a ) + J^T^E^-XiT^Fiit-^dT^ Xi(x) 



(58) 

with 

c i = (u ,X i ), Fi{t) = (F,Xi). 



Using the same technique as in the homogeneous case, it can be shown that 
(|58p is the generalized solution to the problem that can be interpreted as 
the solution in the classical sense under certain additional conditions. 

Let us now consider asymptotics of the solution (|58p in the case of the 
homogeneous initial condition Uo(x) = 0, < x < I and for the source 
function in the form F(x,t) = ip(t) X k (x), X k (x), k = 1,2,... being an 
eigenfunction corresponding to an eigenvalue X k of the eigenvalue problem 
(|2ip - (|22p . The Fourier coefficients of the source function F can be easily 
evaluated: 

Fi(t) = (F, X t ) = (mx k (x), X { ) = rP(t)S ik , 
5ik being the Kronecker delta symbol. The solution (j58j) takes then the form 

u(x, t) = (J T a - l E a ^{-\ k T a )^{t - t) dr^j X k {x), (59) 

so we can analyze its behavior for a given function t/j. To illustrate the 
procedure, let us consider the case if>{t) = tP- l E a ^{-\t a ), a,/3,A > 0. 
Then for A ^ \ k the integral in (|59p can be evaluated in a closed form (see 
[24] ) as the Laplace convolution of two Mittag-Leffler functions: 

u(x,t) = XE ^M-^ a )-X k E a ^(-X k t^ ta+ ^ 1Xk{x ^ (6Q) 

A — Afc 

To extend this result for the case A = X k , we have to differentiate the 
function /(A) = XE a>a+ p(— Xt a ) with respect to A and evaluate the result 
in the point A = A^. The derivative of the Mittag-Leffler function can be 
determined using the term by term differentiation that leads to the formula 

(x^O) 

Using this result, we can explicitly write down the solution u in the case 
A = Afe: 

u(x,t) = i( J E aiQ+ ^_ 1 (-A fe r) + (l-/3)S Qia+/3 (-A fc r))^- 1 X fc (x). 

(61) 

Even if the formulas (j60|) and ()61 j) for the solution u look different in the 
cases A / At and A = A&, the asymptotics of solution as t — > +oo is the 
same in these two cases. To show this fact, the well-known formula for 
asymptotics of the Mittag-Leffler function (see e.g. [23] ) 

Eaj}(-x) = V }~ X) * + OCX" 1 ""), x -> +oo, p = 1, 2, . . . . 



is used. Applying this formula to the Mittag-Leffler functions in (|60p and 
([ST]) , we get the following result: 

Theorem 7 The solution to the initial-boundary-value problem V35 \ )-[31 ty 
for the one- dimensional time- fractional diffusion equation A30\) with the ho- 
mogeneous initial and boundary conditions and with the source function in 
the form 

F{x,t)=i>(t)X k {x), iP(t)=tP- 1 E a:(3 (-\t a ), a,/3,A>0 

is given by A60\) in the case A ^ A& and by \61\) in the case A = A&. The 
asymptotics of the solution u(x,t) = <f){t)Xk{x) as t —> +00 is the same in 
both cases: 




0{tP- a ~ l ), t -> +00, ifa^fi, 
O^""" 1 ), if a = 0. 



5 Conclusions and open problems 

In the paper, some initial-boundary-value problems with the Dirichlet 
boundary conditions for the time-fractional diffusion equation were consid- 
ered. Of course, the same method can be applied for the initial-boundary- 
value problems with the Neumann, Robin, or mixed boundary conditions. 

The maximum principle enables us to obtain information regarding so- 
lutions of differential equations and a priori estimates for them without any 
explicit knowledge of the form of the solutions themselves, and thus is a 
valuable tool in scientific research. In the paper, a maximum principle for 
the time-fractional diffusion equation ([1]) was discussed and applied for prov- 
ing uniqueness of solution of the initial-boundary- value problem ([4])- ([5]) for 
equation ([1]). Of course, following the lines of applications of the maximum 
principle for the parabolic and elliptic PDEs (see e.g. the recent book [27]). 
a lot of other properties of solutions to the time-fractional PDEs can be 
established. In particular, the maximum principle can be applied for some 
classes of the non-linear equations of the fractional order, too. 

Another important and interesting problem that is still waiting for its 
solution would be to try to extend the maximum principle to the space- and 
time-space- fractional PDEs. These equations are very actively employed 
nowadays in modeling of relevant complex phenomena like anomalous diffu- 
sion in inhomogeneous and porous mediums, Levy processes and Levy flights 
and the so called fractional kinetics. Similar to the case of time- fractional 
equations, several different definitions of the space- fractional derivatives are 



used in these equations. From the physical viewpoint, an appropriate maxi- 
mum principle is expected to be fulfilled for the corresponding models. Clear 
understanding regarding what definition of the space-fractional derivatives 
enables a maximum principle would essentially help in attempts towards 
modeling of the real phenomena with the "right" space- and time-space- 
fractional PDEs. 

As to the existence of solution to initial-boundary-value problems for the 
time-fractional diffusion equation, the method presented in the last section 
can be extended for the n-dimensional problems, too. 
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